PROCEEDINGS 

0F SCIENCE 



MCRG study of 8 and 12 fundamental flavors 



Gregory Petropoulos , Anqi Cheng, Anna Hasenfratz and David Schaich 

Department of Physics, University of Colorado, Boulder, CO 80309 
Email: gregory .petropoulosScolorado . edu 

o 

Q We study the renormalization group properties of SU(3) gauge theories with Nf = 8 and 12 nearly- 

y— I massless fermions, using Monte Carlo Renormalization Group (MCRG) two-lattice matching 

i | techniques to predict bare step-scaling functions Sb- Traditional MCRG two-lattice matching re- 

■i—i 

C$ quires that the renormalization scheme be optimized for each bare lattice coupling, so that sy 

I is a composite of many different discrete B functions. We propose an improved procedure that 

Ph 

uses the Wilson flow to eliminate the need for this optimization of the RG blocking transforma- 
tion. While our 12-flavor results indicate an infrared fixed point, si, for Nf — 8 is significantly 
different from zero until strong-coupling lattice artifacts obstruct two-lattice matching. Although 
both procedures produce qualitatively similar bare step-scaling functions, the new s/, obtained 
by combining the Wilson flow with MCRG two-lattice matching have the distinct advantage of 
corresponding to unique discrete j3 functions. 
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1. Introduction and overview of methods & results 

In recent years, many groups have initiated lattice investigations of strongly-coupled gauge- 
fermion systems beyond QCD. While the ultimate goal of these efforts is to explore potential new 
physics beyond the standard model, an essential step is to improve our theoretical understanding 
of the basic properties of these non-perturbative systems. Here we study the renormalization group 
properties of SU(3) gauge theories with Nf = 8 and 12 nearly-massless fermions in the fundamental 
representation, through the Monte Carlo Renormalization Group (MCRG) two-lattice matching 
technique. This is one of several complementary analyses we are currently carrying out, two more 
of which (investigating Dirac eigenmode scaling and finite-temperature transitions) are discussed 
in other contributions to these proceedings [1,2]. Recent references on SU(3) gauge theories with 
Nf = 8 and 12 include [3, 4, 5, 6, 7]; earlier works are reviewed in Ref. [8]. 

In Refs. [9, 10], one of us studied MCRG two-lattice matching for the 12-flavor system with 
nHYP-smeared staggered actions very similar to those we use here. Our gauge action includes both 
fundamental and adjoint plaquette terms, with coefficients related by /3a = — 0.25/Jf. The negative 
adjoint plaquette term lets us avoid a well-known spurious ultraviolet fixed point caused by lattice 
artifacts, and implies jSf = 12/g 2 at the perturbative level. In our fermion action, we use nHYP 
smearing with parameters (0.5,0.5,0.4), instead of the (0.75,0.6,0.3) used by Refs. [9, 10]. By 
changing the nHYP-smearing parameters in this way, we can access stronger couplings without 
encountering numerical problems. At such strong couplings, for both Nf = 8 and Nf = 12 we 
observe a lattice phase in which the single-site shift symmetry ("S 4 ") of the staggered action is 
spontaneously broken [11, 2]. 1 In this work we only investigate couplings weak enough to 

avoid the ^lattice phase. 

In the next section, we review how the MCRG two-lattice matching technique determines 
the step-scaling function Sb in the bare parameter space. Although working entirely with bare 
parameters would be disadvantageous if our aim were to produce renormalized phenomenological 
predictions for comparison with experiment, our current explorations of the phase structures of 
the 8- and 12-flavor systems benefit from this fully non-perturbative RG approach, especially for 
relatively strong couplings. In Section 3. 1 we present our results from the traditional MCRG two- 
lattice matching technique. While our 8-flavor Sj, is significantly different from zero, for Nf = 12 
we observe Sb < for jS^ < 8, indicating an infrared fixed point (IRFP). 

We emphasize that while the existence of an IRFP is physical (scheme-independent), the cou- 
pling at which it is located depends on the choice of renormalization scheme. A limitation of 
traditional MCRG two-lattice matching is the need to optimize the RG blocking transformation 
separately for each lattice coupling j8f- As we explain below, this optimization forces us to probe 
a different renormalization scheme for each so that the bare step-scaling function we obtain is 
a composite of many different discrete /3 functions. 

To address this issue, in Section 2.2 we propose a new, improved procedure that predicts a 
bare step-scaling function corresponding to a unique j8 function. This improved procedure applies 
the Wilson flow [12, 13] to the lattice system before performing the RG blocking transformation. 
Because the Wilson flow moves the system in the infinite-dimensional space of lattice-action terms 

'Ref. [6] recently interpreted the ^lattice phase in terms of relevant next-to-nearest neighbor interactions. 
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without changing the lattice scale, we can use it to approach the renormalized trajectory corre- 
sponding to a fixed RG blocking transformation. By optimizing the flow time tf at each coupling, 
all with the same renormalization scheme, we can carry out the two-lattice matching without a 
need for further optimization. We present some promising but preliminary results of this approach 
in Section 3.2. 

2. Two-lattice matching procedures and the need for optimization 

Two-lattice matching is most easily described in the context of confining systems, where it 
locates pairs of couplings (j3f,j3£) for which lattice correlation lengths obey B, (j8f) = 2^ (j3£). We 
proceed by repeatedly applying RG blocking transformations (with scale factor s = 2) to lattices 
of volume 24 3 x48, 12 3 x24 and 6 3 x 12. 2 Under RG blocking on the m = critical surface, the 
system flows toward the renormalized trajectory in irrelevant directions, and along it in the relevant 
direction. By blocking the larger lattices (with /3f) ii}, times and the smaller lattices (with j3£) only 
rib — I times, we obtain blocked systems with the same lattice volume. If these blocked systems 
have both flowed to the same point on the renormalized trajectory, then we can conclude that 
% (Pf) = 2| (fi' F ) on the unblocked systems, as desired. 

We determine whether the blocked systems have flowed to the same point on the renormalized 
trajectory by matching several short-range gauge observables: the plaquette, all three six-link loops, 
and two planar eight-link loops. For a given j8f , each observable may predict a different AfHp = 
Pf — P'f- The spread in these results is a systematic error that dominates our uncertainties. 

In an IR-conformal system, the gauge coupling that is relevant at the perturbative gaussian 
FP becomes irrelevant at the IRFR The renormalized trajectory connects these two fixed points. 
When RG flows approach this renormalized trajectory, the two-lattice matching can be performed 
and interpreted the same way as in confining systems. In this region the gauge coupling flows to 
stronger couplings, Aj8f > corresponding to a negative RG j3 function. The situation is less clear 
at stronger couplings where we might naively expect backward flow. If there is no ultraviolet FP 
in this region to drive the RG flow along a renormalized trajectory, the two-lattice matching might 
become meaningless. This issue affects every method that attempts to determine the flow of the 
gauge coupling in IR-conformal systems at strong coupling. In all published studies that report an 
IRFP, backward flow has only been observed in a very limited range of couplings in the immediate 
vicinity of the IRFP (cf. Ref. [8]). 

Since we can block our lattices only a few times, we must optimize the two-lattice matching by 
requiring that consecutive RG blocking steps yield the same A/3f . We identify the optimized A/3f 
with the bare step-scaling function Sb- In the following subsections, we describe two different ways 
to perform this optimization. The traditional technique optimizes the RG blocking transformation 
(renormalization scheme). The new method we propose in Section 2.2 instead applies the Wilson 
flow to the lattice system prior to RG blocking, and optimizes the flow time tf. 

2.1 Traditional MCRG 

As in Refs. [9, 10], we use RG blocking transformations that include two sequential HYP 
smearings with parameters (a, 0.2, 0.2), and optimize a as shown in the left panel of Fig. 1. Qual- 

2 We are currently generating larger lattices up to 32 3 x 64, which will permit additional consistency checks. 
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itatively, this optimization finds the renormalization scheme for which the renormalized trajectory 
passes as close as possible to the lattice system with coupling fip. Without optimization, residual 
flows in irrelevant directions can distort the results: this is the reason Afip changes with a in Fig. 1, 
and also explains why increasing the number of blocking steps reduces this a -dependence. 

The downside of optimizing the RG blocking transformation in this manner is that we have to 
use a different renormalization scheme for each jip . As a result, the bare step-scaling function we 
obtain is a composite of many different discrete j8 functions. 



Figure 1: Examples of two-lattice matching optimization for 12-flavor systems. Left: Optimization of the 
HYP-smearing parameter a in the RG blocking transformation, for j3/r = 5.0. Right: Optimization of the 
Wilson flow time tf with fixed a = 0.5, for /3/r = 4.5. In both cases, the uncertainties on the data points are 
dominated by averaging over the different observables as described in the text. 

2.2 Wilson-flowed MCRG 

As an alternative to optimizing the RG blocking transformation, and thus changing the renor- 
malization scheme at each coupling f}p, here we propose to use the Wilson flow to move the lattice 
system as close as possible to the renormalized trajectory of a fixed renormalization scheme. 

The Wilson flow is a continuous smearing transformation [12] that can be related to the MS 
running coupling in perturbation theory [13]. Refs. [14, 15] recently used the Wilson flow to 
compute a renormalized step-scaling function in a way similar to Schrodinger functional methods. 
While this approach appears very promising, it is based on perturbative relations that are only fully 
reliable at weak coupling. Here we do not use this perturbative connection, instead applying the 
Wilson flow as a continuous smearing that removes UV fluctuations. The Wilson flow moves the 
system along a surface of constant lattice scale in the infinite-dimensional action-space; it is not a 
renormalization group transformation and does not change the IR properties of the system. 

Our goal is to use a one-parameter Wilson flow transformation to move the lattice system 
as close as possible to the renormalized trajectory of our fixed RG blocking transformation. We 
proceed by carrying out two-lattice matching after applying the Wilson flow for a flow time tf on 
all lattice volumes. (The Wilson flow is run only on the unblocked lattices, not in between RG 
blocking steps.) As above, since we can block our lattices only a few times, we must optimize tf 
by requiring that consecutive RG blocking steps yield the same Afip , as shown in the right panel of 
Fig. 1 . As for traditional MCRG, increasing the number of blocking steps reduces the dependence 
on the optimization parameter; in the limit rib —> °°> °ur results would be independent of tf. 

With Wilson- flowed MCRG we can efficiently determine bare step-scaling functions that cor- 
respond to unique RG /3 functions. This new capability opens up interesting directions for future 
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studies. By comparing different /3 functions around the perturbative gaussian FP, we can study 
scaling violations in the lattice system. In IR-conformal systems, we can investigate the scheme- 
dependence of the [5 function near the IRFP, an issue explored in perturbation theory by Ref. [16]. 

3. Results 

3.1 Traditional MCRG 

Our results for the bare step-scaling function Sb from traditional MCRG two-lattice matching 
are shown in Fig. 2. On the largest 24 3 x48 lattices that we use in this current study, we work with 
fermion masses m = 0.0025 to stay near the m = critical surface. Under RG blocking with scale 
factor s, the fermion mass changes as s l+r '" where y m is the mass anomalous dimension. Therefore 
we use m = 0.01 on 12 3 x24 and m = 0.02 on 6 3 x 12 lattices. We have explicitly checked that these 
masses are small enough to introduce only negligible finite-mass effects, by generating lattices with 
m = for some points and obtaining indistinguishable results. 

While our 8-flavor results for si, are significantly different from zero for all couplings we can 
explore, for Nf = 12 we find < for j8f < 8, indicating an IRFR Recall that our optimization 
of the RG blocking transformation means that we use a different renormalization scheme for each 
coupling f}p, so these bare step scaling functions are composites of several different discrete /3 
functions. For example, with Nf = 12 at 5 < ftp < 6, our optimization selects renormalization 
schemes with the fixed point near j8f , so that Sb is roughly consistent with zero over an extended 
range. Both our Nf = 8 and 12 simulations encounter the ^lattice phase at strong coupling, where 
we cannot perform matching. As in Refs. [9, 10], we do not explore weak enough couplings to 
recover the two-loop perturbative predictions Sf, « 0.6 (0.3) for Nf = 8 (12). 

As mentioned above, the error bars shown in Fig. 2 are dominated by the spread in results 
from matching four-, six- and eight-link loops. Each of these observables can predict a different 
optimal a, and for fixed a each can predict a different Aj3f. Preliminary results presented at the 
conference determined uncertainties from the full spread of optimal a predicted by the different 
observables. Here, instead, we average A/3f for the different observables at fixed a, and use these 
combined data to optimize a and find the associated uncertainties. 



N f = S 



Two-loop 



5.5 6 

Pf 




Figure 2: Results for the bare step-scaling function st, from ttaditional MCRG two-lattice matching with 
24 3 x48, 12 3 x24 and 6 3 x 12 lattice volumes, for Nf — 8 (left) and Nf = 12 (right). The blue dashed lines 
are perturbative predictions for asymptotically weak coupling. 
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3.2 Wilson-flowed MCRG 

Fig. 3 presents our results for the bare step-scaling function Sb from Wilson-flowed MCRG 
two-lattice matching. We continue using two sequential HYP smearings in our RG blocking trans- 
formation, but now fix the smearing parameters to (0.5,0.2,0.2). We again use24 3 x48 and 12 3 x24 
lattices with fermion masses m = 0.0025 and m = 0.01, respectively, and determine uncertainties 
in the same way as described in the previous subsection. In this preliminary study we don't yet 
employ the volume-corrected optimization discussed in Ref . [9] ; for our lattice volumes, Ref . [9] 
found that neglecting this finite-volume correction introduces only a small additional error. Our 
final results will use the appropriate optimization, and will also present further consistency checks 
from larger lattices up to 32 3 x 64. 
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Figure 3: Preliminary results for the bare step-scaling function Sj, from Wilson-flowed MCRG two-lattice 
matching with 24 3 x 48 and 12 3 x 24 lattice volumes, for Nf = 8 (left) and Nf = 12 (right) with fixed HYP- 
smearing parameters (0.5,0.2,0.2). The blue dashed lines are perturbative predictions for asymptotically 
weak coupling. 



While our results in Fig. 3 from combining the Wilson flow with MCRG two-lattice matching 
are qualitatively similar to the results of the traditional approach in Fig. 2, we can now identify the 
bare step-scaling function st with a unique discrete /3 function. In the renormalization scheme de- 
fined by our RG blocking transformation with HYP-smearing parameters (0.5,0.2,0.2), we find a 
12-flavor IRFP at 5 < ftp < 6. Although fip is scheme-dependent, the existence of this IRFP is phys- 
ical. We are currently exploring other choices of renormalization schemes, to non-perturbatively 
investigate the scheme-dependence of the j3 function near the IRFP [16]. At weaker couplings, we 
will also attempt to use similar explorations to study scaling violations in our lattice systems. 



4. Conclusions 



We have proposed a new, improved MCRG two-lattice matching procedure that uses the Wil- 
son flow to eliminate the need for optimization of the RG blocking transformation. Both traditional 
MCRG and Wilson-flowed MCRG produce bare step scaling functions Sb that indicate an infrared 
fixed point for SU(3) gauge theory with Nf = 12 fundamental fermions, while Sb for Nf = 8 is 
significantly different from zero in the accessible range of lattice couplings. The results obtained 
by combining the Wilson flow with two-lattice matching correspond to a unique /3 function, unlike 
Sb from traditional MCRG, which is a composite of many different discrete /3 functions. 
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